Dynamical polarization, screening, and plasmons in gapped graphene 
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The one-loop polarization function of graphene has been calculated at zero temperature for 
arbitrary wavevector, frequency, chemical potential (doping), and band gap. The result is expressed 
in terms of elementary functions and is used to find the dispersion of the plasmon mode and the 
static screening within the random phase approximation. At long wavelengths the usual square root 
behaviour of plasmon spectra for two-dimensional (2D) systems is obtained. The presence of a small 
(compared to a chemical potential) gap leads to the appearance of a new undamped plasmon mode. 
At greater values of the gap this mode merges with the long-wavelength one, and vanishes when 
the Fermi level enters the gap. The screening of charged impurities at large distances differs from 
that in gapless graphene by slower decay of Friedel oscillations (1/r 2 instead of l/r 3 ), similarly to 
conventional 2D systems. 

PACS numbers: 73.21.-b, 73.61.Wp, 81.05.Uw 



I. INTRODUCTION 

Graphene, a free-standing flat single layer of carbon 
atoms, has attracted considerable attention since its 
experimental discovery [l[ due to its unique electronic 
properties. The conduction and the valence bands 
in graphene touch each other at zero energy at two 
inequivalent points, K and K', in the momentum 
space. Close to those points charge carriers have a 
linear isotropic energy spectrum £7±(k) = ±7iw F |k| 
and behave like relativistic massless Dirac particles 
in 2 + 1 dimensions, with the Fermi velocity vp ~ 
10 6 m/s playing the role of the speed of light Q. This 
remarkable feature makes graphene different from all 
other known two-dimensional systems with their usual 
quadratic dispersion law of charged quasiparticles. For 
the undoped graphene at zero temperature the valence 
band is completely filled ('Dirac sea') and the conduction 
band is empty, which corresponds to a zero value of the 
chemical potential (Fermi energy) \x. The latter can 
be made non-zero and regulated experimentally by the 
applied external electric field [l[ 

One of the most important fundamental quantities 
for understanding physical properties of graphene is the 
dynamical polarization function II(cj, k) which describes 
the screening of the Coulomb potential due to many- 
body effects and determines the collective excitation 
modes. This function for graphene at finite chemical 
potential and u> = has been calculated in the one-loop 
approximation in Later, in Q and Q it was obtained 
for the arbitrary frequency, and the results were used, in 
particular, to find the dispersion of the plasmon mode 
within the random phase approximation (RPA). 

The present paper deals with the more general 
case, when a finite gap 2A between the conduction 
and the valence band exists in graphene, and the 
charge carriers behave analogously to relativistic massive 
Dirac particles, with the dispersion relation i?±(k) = 

A 2 . Such a gap is observed in experiments 



graphite [8| substrate (with the band gap 2A of about 
0.26 eV or 10 meV, respectively), and is commonly 
attributed to the sublattice symmetry breaking due 
to the commensurate perturbation from the substrate. 
Moreover, even in the absence of any substrate, there 
should still remain a small gap (2A so ~ 10~ 3 meV) 
originating from the spin-orbit interaction Q. The 
gap can also be generated dynamically by applying an 
external magnetic field [3, [l(|. The dispersion of a 
plasmon mode for graphene with the finite spin-orbit 
fermionic gap has been obtained numerically in [ll( . Our 
aim is to explore in more detail the influence of the gap on 
the plasmon dispersion and the static charge screening. 
In this work we restrict ourselves to the case of zero 
temperature and calculate analytically the dynamical 
polarization function for gapped graphene. Next we 
find the screened Coulomb potential at large distances 
and the spectrum of plasmons for different values of the 
dimensionless parameter A//x by solving the dispersion 
equation with the obtained polarization function. 



II. MODEL 

The low-energy Hamiltonian of graphene can be 
written in a Dirac- like form 12], 

Ho = -ihv F ^2 *1( T 3 ® <ridi +T (g) cr 2 9 2 )* s 



s=± 



-ifivp ^si^di + 7 2 9 2 )* s 



(1) 



s=± 



± 



lfi 2 v 2 F k 2 



where * s = {Vka^kb^k'A^k'bY is a four - 
component wavefunction of charged quasiparticles, the 
components of which refer to the A and B sublattices and 
to two inequivalent valleys in momentum space (around 
K and K' points). Pauli matrices Oi and t; act on 
sublattice and valley components, respectively (<7q, tq 
are unit matrices), s is a spin index, gamma matrices are 
defined as 7 = r ® 03, 7 1 = zt 3 ® er 2 , 7 2 = — ir ® a\, 



with graphene placed on the silicon carbide 7] or and ^> s = ^^7° is the Dirac conjugated spinor 
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Interaction with a substrate of a certain type may 
break the sublattice symmetry. In this case the 
Hamiltonian |T]) acquires an additional term 

Ha = A J2 *l r o ® = A ^ * s * s , (2) 

which is equivalent to the usual Dirac mass term (in the 
current representation of 7-matrices), and the energy gap 
2A opens between the valence and the conduction band. 
Another gap is driven by a spin-orbit interaction which 
adds the term 

'kCso — A so 

s=± 

to the Hamiltonian. 

III. POLARIZATION FUNCTION 



In the absence of a band gap the expression for the 
polarization function (|A.20[) (|A.21[> reduces to the one 
calculated in [f| and Q , and can be written in the 
following way 

n( „, /:) - 2e2iV /M e*y 2 

e$h vp Ae^fiyJ vpk 2 — lu 2 

x - 2,1 + hu r s 2y + hu y | iaiccozh ( ^ + ^ \ 

fiVFk V ^ hvpk J V hvpk J 

2/j, — huj j /2/i — Tilu\ 2 ^/2/i — 

hvpk V ^ fivpk ) \ hvpk J 

(7) 

with the prescription lu — > lu + iO. This expression is an 
analytic function of lu without singularities in the whole 
upper complex half-plane. Hence it coincides also with 
the result obtained in 

.13] 

when calculated on the positive 

imaginary half-axis. 



The screening of an instantaneous bare Coulomb 
interaction 



U (t,r) 



e 2 S(t) 



duj f d 2 k exp(ikr — iuit) 



2vr / 2ir 



(4) 



due to many-body effects is determined by the retarded 
polarization function U(lu, k) or, equivalently, by the 
dielectric function e(u>, k) = 1 + n(w, k)/k, resulting in 



U(t,r) 



e 
£0 



2^ 



d 2 k exp(ikr — iujt) 

~2ir k + n(u;,k) 



(5) 



where e$ is an effective background dielectric constant, 
the average of that of the substrate and that of the 
vacuum. Note that H(uj, k) is defined here in the same 
way as in [4 , which differs from the polarization function 
used in 

MM 13 by a factor of -2ire 2 /e . 

We consider clean graphene, i.e. the effect of scattering 
of charge carriers on impurities is not taken into account. 
The polarization function is proportional (with a factor of 
2tt/£q) to the time component of the vacuum polarization 
tensor in a three-dimensional quantum electrodynamics 
(QED3), where vf is the 'speed of light' and A/vjr is a 
fcrmion mass. We calculate IL(u>, k) analytically in the 
one-loop approximation. The exact results are given in 
the appendix by the equation (|A.16j) for the case < A 
and by the equations (|A.20|H|A.21j) for fi > A. In what 
follows it is assumed that fi > and lu > 0, taking into 
account that the polarization function depends only on 
the absolute value of \x and U(—lu, k) = U(oj, k)*. 

At small energies and momenta TivfU <C Tilu <C /x and 
fi > A it is easy to obtain from (|A.20[) and (|A.21|) 



n(w,fc) 



e 2 N f k 2 fi ( _ A^ 
M 2 



eoh lu 



2 , ,2 



(0) 



Here Nf is the number of fermions in QED3 (in the case 
of graphene Nf = 2 due to the spin degeneracy). 



IV. STATIC SCREENING 

In the RPA approximation, the static screening of the 
Coulomb potential (j)o(r) = Ze/r of an impurity with 
charge density no(r) = ZeS(r) is determined by the static 
(lu = 0) one-loop polarization function. The induced 
charge density in the graphene plane reads 



6n(r) = Ze 



d 2 k 
(2n) 2 



£o£(0,k) 



exp(ikr) 



and the resulting screened potential is 



. . Ze f d 2 k exp(ikr) 
(r) = — 

£0 



2ir fc + n(0,fc) 



(8) 



(9) 



At zero frequency and /j > A we recover from (jA.20[> — 
(jA~2l|) the previous result for the static polarization 
function d 



n(o,fe) 



2e 2 N f n 
eoh 2 Vp 
h 2 v 2 F k 2 



1 - e(k - 2k F ) 



2k 



AfiVFkp 



4 A 2 hvpx/k 2 
■ arctan — 



2/i 



(10) 



At large distances from the charged impurity (kpr 3> 1), 
there are two main contributions to ([8} and Q. The 
first, Thomas-Fermi contribution, which is determined by 
the long- wavelength (k — > 0) behavior of the polarization 
function, 



n(0,fc < 2k F ) = 



2e 2 N f fi 
e h 2 v 2 F 



is the same as in gapless graphene [j| 
e 2 N f fi 



8riTF{r) 



7TEoh 2 vp 



■4>TF(r) = - 



Zh 2 v 2 F 



(11) 



(12) 
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and decays under the same 1/r 3 law as in an 
ordinary non-relativistic two-dimensional electron gas 
(2DEG) [14J. The second contribution is oscillatory 
(Friedel oscillations) and comes from the non-analyticity 
of the polarization function ([TO]) at k = 2k F where 
its derivative is discontinuous (for A ^ 0). These 
oscillations 

k F , . , Ze 3 N f h 2 vlk F A 2 sm(2k F r) 
Sn osc (r) = — 4> osc {r) ~ 



2nn(e h 2 v 2 F k F + e 2 N f /j) r 2 v F = 10 6 m/s. 



( 13 ) 

decay as 1/r 2 as in a 2DEG [l4( and slower than in 
gapless graphene where they obey the 1/r 3 law [f| [l5[ . 
The reason for this difference is that at A = the 
polarization function at k = 2k F has a discontinuity only 
in the second derivative, in contrast to the case of the 
gapped graphene and 2DEG. 

At jjL < A, the static limit of (|A.16[) reads 



n(o,jfe) 



e 2 N f ( h 2 v 2 k 2 -4A 2 
J I 2A ' F 



2e h 2 v 2 F \ " hv F k 
hv F k 



(14) 



x arcsm 



h 2 vlk 2 + 4A 2 , 



(compare with equation (14) in [J]), and the induced 
charge density ([8]) at large distances, which is determined 
by the long- wavelength behaviour of (fl4|) , 



11(0, fc^ 0) 



e 2 N f k 2 
3e A ' 



decays as 1/r 3 [lq ]. 



Sn(r) 



Ze 3 N f 
67r£nAr 3 ' 



The potential at large distances for this case 

1) 



Ze ( e 4 N 2 



(r) ~ e r 1 1 9e 2 A 2 r 2 



(15) 



(16) 



(17) 



in the main order remains screened only by the substrate. 



V. PLASMON DISPERSION 

The plasmon dispersion u> p (k) is obtained within the 
RPA approximation by finding zeros of the dielectric 
function e(u>,k) or, equivalently, by solving the equation 



U{oj p ,k) +k = 0. 



(18) 



Using the approximate expression Q we obtain 
immediately for small to and k 



oj p {k) 



e h 



(19) 



with a s quar e root behaviour usual for two-dimensional 
systems |14| . which differs from the case of gapless 
graphene [1, 0] only by the factor y/l — A 2 j \i 2 . 

To obtain the plasmon dispersion for an arbitrary 
frequency and wavevector, the equation (|18p with the 
polarization function defined at (|A.20[) and (IA.21|) is 
solved numerically. We use the following parameters: 
£o = 2.5 (this value corresponds to the case when a 
graphene sheet is deposed on a SiC>2 substrate) and 



Regions 1A-2A and 2B-4B (Fig. |A~T|) in the (k,uj) 
space, in which the imaginary part of the polarization 
function is different from zero, form the single-particle 
excitation (SPE) continuum, displayed in Fig. Q] as a 
shaded area. In those regions plasmons decay into 
electron-hole pairs (Landau damping) and equation l|18p 
has no real solutions. Then one has to solve the complex 
dispersion equation 



H(u)p — «7, k) + k = , 



(20) 



finding both energy Tuj p {k) and the decay rate 7(fe) of 
plasmons. We solve (f20f numerically, using analytical 
continuation of (|A.20[1 — (|A.21[1 from the real to axis into 
the lower complex half-plane. 

When the damping of the collective mode is weak, the 
approximate equation 

3? e n(w p ,fc) + fc = (21) 

is often solved [H, [ll[ instead of (f2"0"|) and the decay rate, 
which is assumed to be small, is determined from another 
equation, 



3m n(w p , fc) 



{d/duj)Vlen(u],k)\ 



(22) 



The numerical solutions of equations ([20]) (for the 
frequency and the damping rate) and (|21~j) together with 
the approximate analytical solution (|19p are plotted in 
Fig. [U for different values of the gap parameter. 

At zero gap there is one plasmon mode which is 
undamped until it enters the SPE region (Fig. [TJa)) [1, 
6]. When the gap is different from zero but is still 
small (A < 0.22/i), a new undamped plasmon mode 
emerges in the gap which opens between two SPE regions 
(corresponding to intra- and interband transitions of 
electrons) in the vicinity of the line u> — v F k, as displayed 
in Fig.[ljb). The damping part of the plasmon spectrum 
is no longer a continuous extension of the undamped part 
but rather a separate mode (Fig.QJe)). When A becomes 
larger than ~ 0.22/i (this threshold value depends on 
the background dielectric constant, for example in the 
vacuum, when £q = 1, it is approximately equal to 0.29/i), 
the damped part of the spectrum disappears and two 
undamped modes merge (Fig. [He)). With the further 
increasing of A, the plasmon dispersion curve becomes 
shorter (Fig. [Ud)) and vanishes at A = fi. The gap 
dependence of the plasmon spectrum described above is 
also shown in Fig. [51 



1 ■ ' ' ' ■ ' 1 

1 2 3 1 2 3 




hv F k/p. hv F k/p. 

FIG. 1: Panels (a), (b), (c), and (d) show the plasmon spectrum for A = 0, A = 0.2/i, A = 0.3/j,, and A = 0.9/i, respectively, 
and (e) is a magnified fragment of (b). Solid lines represent the solutions of the complex dispersion equation, dotted lines 
are the solutions of the approximate equation (|2ip and dashed curves show the long-wavelength square root dependence. The 
shaded area indicates the SPE regions in which plasmons decay with the damping rate shown in panel (f) for A = (dashed 
line) and A = 0.2/i (solid line). 



Thus wc found that plasmons are absent at A > /i. 
Indeed, the dispersion equation (fT8j) can be satisfied only 
when IIo(w, k) < 0, and it can be shown from (IA.15|) 
that ^teHo(u>, k) > in the region where Ho(ui,k) is real 
(however, for the gapless case this conclusion was shown 
to be an artefact of the RPA approximation [l7|, and 
the ladder corrections to polarization function lead to 
the existence of the undamped plasmon mode even at 
fi = A = 0). 



In those regions where the solutions of equations (|20|) 
and (f2"Tj) are close to each other, the decay rate given 
by (|22[) is small and does not differ significantly from 
7 obtained from (I20[) . When the solutions move away 
from each other, it turns out that the expression (|22p is 
not small, and the approximate equations (f2"Tj) and (l2"2")) 
are not valid. In particular, the lower branch which 
appears for A ^ and lies fully in the intraband 
SPE region (Figs, [ljb)-(e)) is merely an artefact of this 
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FIG. 2: Coordinates (momenta) of the undamped plasmon 
mode edges as a function of the gap value: the end of the 
mode starting at k — (solid line), the beginning (dotted 
line) and the end (dashed line) of the second mode which is 
absent in the gapless case. 
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approximation. 

Our numerical results are in agreement with the 
work [ll|, where the equation (|2Tj) has been solved and 
a graph, similar to Fig. [He) was obtained for the case of 
zero temperature. 



VI. CONCLUSIONS 

We have calculated analytically the one-loop 
polarization function in the presence of a finite gap in the 
spectrum of Dirac quasiparticles. The result is expressed 
in terms of elementary functions and is valid for arbitrary 
wavevector, frequency, doping, and gap. The obtained 
expression is employed to calculate the static screening 
of Coulomb potential and the spectrum of plasmons 
within the RPA approximation. The oscillating part of 
the screened potential decays slower than in the gapless 
graphene, while the non-oscillating part is the same. 
It turns out that collective plasmon modes are present 
only at A < fi, i.e. when the Fermi level lies outside the 
gap. At long wavelengths the plasmon dispersion shows 
the usual square root dependence for 2D systems. At 
non-zero gap, an additional undamped mode is present 
at small values of A//i, and at 0.22/i < A < fj, only one 
undamped plasmon mode remains. 
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APPENDIX: CALCULATION OF n(w,fc) 

To obtain the retarded polarization function we 
start with a finite temperature T. In the one-loop 



FIG. A.l: Regions with different expressions for polarization 
function (here A = 0.86/i). At the borders H(u,k) is smooth 
(dotted line) and has a logarithmic singularity (solid line) or 
discontinuity of its derivative (dashed line). 



approximation it is given by the expression 

d 2 q 



£ o n t^J ( 27r ) (AT) 
x tr[7 5(jw m + iQ n ,k + q)7oS , (ifi n ,q)] , 

where uj m = 2i:mT '/% and Q n = (2n + \~)kT/% are 
Matsubara frequencies, and the propagator of Dirac 
quasiparticles has the form 



S(iVL ni c\) 



(ihQ, n — ju)7° — hvFCff — A 



(A.2) 



Here we are not taking into account the spin-orbit gap, 
considering only the sublattice symmetry breaking mass 
term @. After taking the trace over spinor indices we 
obtain (for convenience we shall omit constants fi and vf 
during intermediate calculations) 



U(iw m ,k) 



2e 2 TN f 



+ OC 



d 2 q 



(iuj m + ifi n - (J,)(iQn - n) + q(k + q) + A 2 
[(iu m + m n - n) 2 - E 2 +q \ [(in„ - fi) 2 - E 2 ] ' 



(A.3) 



with Eq = \J q 2 + A 2 . To perform the summation over 
n we rewrite the integrand as 



E 

A,A'=± 



4 V S q Sk+ q J 



1 



1 



(A.4) 



iSl n — (j, + XEq ifl n + ioj m — fi + X'Ek +c 
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where A, A' are band indices, and then use the formula 
+ °° 1 1 

i.Q„ — // -I- t*i iO._ -I- id)— — ii -4- To 

(A.5) 



E 



ifl n — fx + x% iVl n + io; m — fx + X2 

1 n F {x\) - uf{x2) 



T X\ — X2 — iw r , 



where n F (x) — \e.( x ~^/ T + l] . After summation, 
the analytic continuation from Matsubara frequencies is 
made by the replacement iu> m — ► cu + iO, and the retarded 
polarization function reads 



n <^>=-S/ d2 < ? ( i+aa 



,q(k + q) + A 2 



n F (\Eg) ~ n F (\'Ek +cl ) 
\Eq - A'£ k+q - a; - iO 



E n E\ 



q^k+q 



(A.6) 



The same expression (with A so instead of A) is obtained 
if the sublattice symmetry is unbroken and the spin-orbit 
gap © is taken into account , hence our results will be 
valid for both cases. At zero temperature Fermi functions 
n F (x) turn into step functions (9(/i — x), and (|A.6[) gives 

1%, k) = -x" (w, k) + fc) + *) , (A.7) 

where 

x | 1 ± 



e 2 A^ 



27T£ 

q(k + q) + A 2 



d 2 qO(D 2 - A 2 - (j 2 



-Eq-E-k+q 
1 



(A. 



1 



w + £; q =F #k+ q + iO uj- Eq± £ k4 



Here the upper and the lower signs correspond 
to intraband and interband electron-hole transitions, 
respectively, and the parameter D defines the integration 
limits. We proceed with further calculations in a similar 
way to those shown in [j| for massless Dirac quasiparticles 
(A = 0). In the case /i < A (when the Fermi level lies in 
the band gap) only the first term in (|A.7[) survives. Thus 
we can split our polarization function into two parts 

U{uj,k) = U (lu, k) + 9(fi - A) Ui(uj, k) , (A.9) 

where IIo(a>, k) = — X^o(w, k) is the polarization function 
at chemical potential /i < A, which does not depend on 
the value of fj,, and 



IIi(u, k) = xt(w, k ) + Xu (w, k) 



(A.10) 



is a correction to the polarization function in the case 
when the chemical potential is greater than the gap 
parameter A. 



At first we consider the imaginary part of (IA.8|) , 



3mx^(w,fc) 



<7 



2e 



9(D — A) J dqq J dip 



x 1 ± 



q(k + q) + A 5 



(ATI) 



EqEk + q 

5{gj + £ q T £k+ q ) - S(u -Eq± £ k+q ) 



Arguments of these 5-functions determine the SPE 
regions in the (k, oj) space. The integration over ip yields 



^rnx D (uJ,k) 



e 2 N 



£0 



6{D - A) J dE 



<3=±A 

W(±l) T S(-^(±E T w) 
(2E + I3u) 2 - k 2 
e((k 2 - oj 2 ) [(2E + M 2 - k 2 ] - 4fc 2 A 2 ) 
(fc 2 - oj 2 ) [{2E + f3u)) 2 - k 2 ] - 4fc 2 A 2 



(A.12) 



Taking the lower sign and setting D = oo in the 
above expression and integrating over E we obtain the 
imaginary part of polarization at [i < A, 



3mn (w,fc) = — Xoo (w, fc) 



4A 2 ) 1 



ne 2 N f k 2 
4e v / w 2 - 

4A 2 
cj 2 - fc 2 



(A.13) 



We find the real part of the polarization function at 
fi < A from the Kramers-Kronig relation, using its 
already known imaginary part (|A.13[) . instead of direct 
calculation of dt,ex^o(oJ, k) from equation (|A.8j) 



oo 

3?eIIo(a;,fc) = — v. p. / dio 

7T J 



,3rn IIo(w', A;) 



OJ' — OJ 



e 2 N f k 2 f , ,fl(cy 2 -£: 2 -4A 2A 
■ v.p. / doj 



4e 



x 1 + 



(oj> ~ oj)Voj' 2 ~ k 2 



(A.14) 



4A 2 



T^ 2 sguuJ 



After integration we have 

e 2 N f k 2 f A 



^eU (oj,k) = 

fc 2 - to 2 ~ 4A 2 r 



£q 1 fc 2 — oj 2 



4 |fc2 _ w 2|3/2 



6*(fc — oj) arccos 



-6{u)-k)hx 



(2A + Vm 2 - fc 2 ) 2 
|w 2 -fc 2 -4A 2 | 



k 2 - oj 2 - 4A 2 
oj 2 - k 2 - 4A 2 



(A.15) 
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We can combine (IA.13[) and (|A.15|) . obtaining the 
retarded version of the time component of the well-known 
expression for the polarization tensor in QED3 (l8l | 



Uo{u,k) 



e 2 N f k 2 



2e Q h(v 2 F k 2 -uu 2 ) 



2 A v 2 F k 2 -uu 2 -4A 2 /h 2 



x arcsm , 



Vpk 2 — LU 2 



v 2 F k 2 -lu 2 + 4A 2 /h 2 J ' 



\Jv 2 F k 2 — oj 2 

(A.16) 



where the prescription lu — > lu + iO is used (compare also 
with equation (6) in fl6(| ) . 

Using the expressions (jA.lOj) and (|A. 12|) we can write 
the correction to the imaginary part of the polarization 
function as 



3m 111 (w, k) 



e 2 N f 



E 



SO — ± JA 



dE 



sgn(u + pE) \(2E + Plu) 2 - f 
9((k 2 - Lu 2 ) [(2E + Puj) 2 - k 2 ] - 4k 2 A 2 ' 
(k 2 - lu 2 ) [{2E + Plu) 2 - k 2 ] - 4fc 2 A 2 



(A.17) 



The real part of Hx(w, k) is obtained from (IA.10|) 
and (|A.8)) . and after some algebraic simplifications it 
reads 



3?eni(w,fc) = 



,2 Nf fV^^ 2 dqq 



tt£o Jo 
(lu + 2£ q )£ q + qk 



E, 



v.p 



q JO 

+ (lu — > -lu) 



(A.18) 



k 2 -Lu(LU + 2Eq) + 2qk 
Integrating this expression over the angle we find 
e 2 N f 



dE 



sgn(fc 2 - lu(2E + uu)) 



^0 UA 

e((uu 2 - k 2 ) [(2E + lu) 2 - k 2 ] + 4k 2 A 2 ) 



k 2 )[(2E- 



k 2 ] + 4fc 2 A 2 



(A.19) 



x [(2E + uu) 2 - k 2 ] +(lu^ -uj) 



2( M -A) 



Performing the last integration in (|A.17|) and (|A.19[) and 
adding the results to (|A. 13|) and (|A. 15|) . respectively, we 
arrive at the final expression for the polarization function 
at fi > A. Restoring constants fi and vp and introducing 
the following notations 



f(u,k) 



2 N f k 2 



4e h^J\v 2 F k 2 



xo = 4 / 1 



4A 2 



h 2 (v 2 F k 2 



LU 2 ) 



G<(x) — x\j Xq — x 2 — (2 — Xq) arccos(a;/xo) , 
G > (x) — x\J x 2 — Xq — (2 — x 2 ) arccosh(a;/a;o) , 
Gq(x) — x\J x 2 — Xq — (2 — Xq) arcsinh/a;/ \J — x 2 ^ 



we can write the results in the next form 
2e 2 N f fj, 



3ieU(Lu,k) 

r 0, 



e h 2 v F 



f(tu,k) 



x < 



G< 



G, 



2/1 — flbJ 

Tivpk 
2/i + Tioj 



V Tivpk J \ Tivwk J 



hvFk 

' 2/i — flLU 



fLVFk 
2jJL + TlLU'> 



t Zfj, - noj \ /Zfi + tiuu\ 

\ fivpk ) \ hvpk ) 



Tivpk 



hvpk 
2/i — Tvjj^ 



//fi + noj\ /Z(j,-hLu\ 
V Tivvk ) \ fivwk J 



hvpk 
2/i + Hlu 
hvpk 



^, '2fi + htu\ g fhu — 2(i 



G~; 



Gt 



hvpk / V hvpk 

/hcu - 2/i\ g /2fj. + hu)\ 
V TlVFk J \ hvpk J 



hvpk 
2/i + Hlu 
hvpk 



G 



•( 



hvpk 
2/i — 

TlVFk 



$smll(w, k) — f(u>, k) 



,2/i + 7kj\ /2fi — hLU 

G>[ ^~^~ ] ' G> {~h^k~ 



x < 



0, 
0. 

0. 



hvpk I 
/ 2/i + ?iw \ 
V Tivpk J 



G. 



/2/i — ftcjx 
V hvpk J 



n(2-x 2 ) 
tt(2-x 2 ) 
10, 



1A 
2A 

3A 

4A 

IB 

2B 

3B 

4B 

5B 



(A.20) 



1A 



2A 
3A 

4A (A.21) 
IB 

2B 

3B 

4B 
5B 



with the following regions in the (k,u>) space (shown in 

Fig.EU) 



1A: Hlu < /i — \/h v F (k — kp) 2 + A 2 



2A: ±/i=F \/h 2 v 2 F (k- kp) 2 + A 2 < Hlu 



H 2 v 2 F (k + kp) 2 + A 2 



< -ft- 



3A : hLu < -fj, + \/h 2 v 2 F (k - kp) 2 + A 2 , 



4A : -fi+ Jh 2 v 2 F (k + k F ) 2 + A 2 < Hlu < hv F k . 



IB: k<2k F , 



H 2 v 2 F k 2 + 4A 2 < Hlu 



<f i+Jn 2 v 2 F (k-k F ) 2 + A 2 , 



8 



2B 



n + ^Jh 2 v 2 F (k 



A 2 < huj 



<[i + x /n 2 v 2 F (k + k F y + A 2 , 



3B : > n+J h 2 v 2 F (k + k F ) 2 + A 2 , 
4B : k>2k F , Jh 2 v 2 F k 2 + 4A 2 < ?kj 



< Ai + \Jfi 2 v 2 F {k - fc F ) 2 + . 



5B : ?iv F fc < Tiw < Jfi 2 v 2 k 2 + 4A 2 , 



where fc F = ^/m 2 A 2 /hv F . 

In the gapless limit Xo becomes equal to 1, regions 4A, 
4B, and 5B vanish, and the expressions (j A.20|) - (| A.21|) 
reduce to ([7]). 



[1] K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, 
Y. Zhang, S.V. Dubonos, I.V. Grigorieva and 
A. A. Firsov, Science 306, 666 (2004); K.S. Novoselov, 
A.K. Geim, S.V. Morozov, D. Jiang, M.I. Katsnelson, 
I.V. Grigorieva, S.V. Dubonos and A. A. Firsov, Nature 
438, 197 (2005); Y. Zhang, Y.-W. Tan, H.L. Stormcr 
and P. Kim, Nature 438, 201 (2005). 

[2] P.R. Wallace, Phys. Rev. 71, 622 (1947). 

[3] D.P. DiVincenzo and E.J. Mele, Phys. Rev. B 29, 1685 
(1984); G.W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 

[4] E.V. Gorbar, V.P. Gusynin, V.A. Miransky and 
LA. Shovkovy, Phys. Rev. B 66, 045108 (2002). 

[5] B. Wunsch, T. Stauber, F. Sols and F. Guinea, 
New J. Phys. 8, 318 (2006). 

[6] E.H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 
(2007). 

[7] SY. Zhou, G.-H. Gweon, A.V. Fedorov, P.N. First, 
W.A. de Heer, D.-H. Lee, F. Guinea, A.H. Castro Neto 

and A. Lanzara, Nat. Mater. 6, 77 (2007). 

[8] G. Li, A. Luican and EY. Andrei. larXiv:0803.4016l 
[9] D. Huertas-Hernando, F. Guinea and A. Brataas, 
Phys. Rev. B 74, 155426 (2006); H. Min, J.E. Hill, 
NA. Sinitsyn, B.R. Sahu, L. Kleinman and 



A.H. MacDonald, Phys. Rev. B 74, 165310 (2006); 

Y. Yao, F. Ye, X.-L. Qi, S.-C. Zhang and Z. Fang, Phys. 

Rev. B 75, 041401 (2007). 
[10] D.V. Khveshchenko, Phys. Rev. Lett. 87, 206401 

(2001); E.V. Gorbar, V.P. Gusynin, V.A. Miransky and 

LA. Shovkovy, Phys. Rev. B 78, 085437 (2008). 
[11] X.F. Wang and T. Chakraborty, Phys. Rev. B 75, 033408 

(2007). 

[12] C.L. Kane and E.J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

[13] Y. Barlas, T. Pereg-Barnea, M. Polini, R. Asgari and 
A.H. MacDonald, Phys. Rev. Lett. 98, 236601 (2007). 

[14] F. Stern, Phys. Rev. Lett. 18, 546 (1967). 

[15] V.V. Cheianov and V.I. Fal'ko, Phys. Rev. Lett. 97, 
226801 (2006). 

[16] V.N. Kotov, V.M. Pereira and B. Uchoa, Phys. Rev. B 

78, 075433 (2008). 
[17] S. Gangadharaiah, A.M. Farid and E.G. Mishchenko, 

Phys. Rev. Lett. 100, 166802 (2008). 
[18] TW. Appelquist, M. Bowick, D. Karabali and 

L.C.R. Wijewardhana, Phys. Rev. D 33, 3704 (1986). 



